Quantum effects on the phase diagram of nuclear-hke systems 
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A Path Integral Monte Carlo method is used to investigate the thermodynamics of nuclear like 
systems. Systems composed of bosons or fermions interracting via a Lennard- Jones potential with 
periodic boundary conditions were simulated and the corresponding phase diagrams are constructed. 
The Path Integral Monte Carlo appears to be a powerfuU tool for investigating quantum effects in 
nuclear multifragmentation phenomena. 
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Highly excited nuclear systems are good "laborato- 
ries" for thermodynamical studies. Due to the van der 
Waals type of the nucleon-nucleon interaction these sys- 
tems are supposed [l[ to exhibit liquid-gas phase tran- 
sitions. The connection is however not straight forward 
due to their small number of constituents, the presence 
of the Coulomb force and quantum effects. The effect of 
the Coulomb interaction on the nuclear liquid-gas phase 
transitions was previously studied (see e.g. Refs. @, 
The effects of other multifragmentation specific features 
such as finite size (for attempts to describe the finite 
size effects of the system within Hartree-Fock theories 
see Refs. 0, and degree of homogeneity were ad- 
dressed in Refs. [6|-l9|. Nuclear multifragmentation of ex- 
cited nuclear systems was treated wi th clus ter-type mod- 

In this kind 



els lyke SMM, MMMC, and MMM 
of models a fragmentation event is approached by various 
size excited nuclear fragments interacting only via their 
Coulomb field. This approach works very well at small 
densities where, the average distance between fragments 
is large enough and fragments don't interact via nuclear 
forces. At larger densities, fragments are nolonger spher- 
ica 1 and the nuclear interaction between them is quite 
important. At such densities, quantum and shape de- 
generacy effects play an important role. Indeed, it was 
proven that a cluster-type statistical multifragmentation 
model cannot accurately describe the dense branch of the 
phase diagram 

Monte Carlo simulations with classical Lennard- Jones 
fluids showed that the accurate treatment of fragments' 
shape degeneracy leads to the restoration of the Guggen- 
heim shape of the phase diagram 0, However, the 
systems considered there were classical. What about 
quantum ones? This question is addressed in the present 
paper. 

In the present work we perform a simulation of a quan- 
tum fluid interacting via a Lennard Jones 6-12 (LJ) fluid. 
Such simulations are tractable with Path Integral Monte 
Carlo (PIMC) methods and have been extensively used 
in sohd-state physics (see eg. (S])- The PIMC method 
is based on the isomorfism that one can acheive between 
the average value (O) of an observable (O) of a N-particle 



quantum system and the average value of the the same 
observable correspnding to a classical N-polymer system. 
This correspondence was firstly pointed by Feynman in 
1972 lj|. Let us briefly present some Path Integral the- 
ory. The average value of a system observable, O, writes: 



(1) 



where Z = ^ - e~^^' is the canonical partition function 
is the ith particle wave-function. In a position-space 
representation, the density matrix writes: 

p(R, R'; (3) = (R|e-^^' |R') = e-^''^'^*^^ (2) 

i 

where R = {ri, r2, Tat}. So: 

{0) = ^ J dRdR'p(R,R';/3) (R|0|R') (3) 

For a free particle in a _D-dimensional box of size L with 
periodic boundary conditions (PBC), one has: 
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Hence: 



p{r, r';P) = -^Yl i^^^^l + *k„(r - r')] 
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where for a particle of mass m, \^n^/2m. The under- 
lying principle of introducing path integrals in imaginary 
time is the product property of the density matrix stating 
that the low temperature density matrix can be expressed 
as a product of high-temperature density matrices: 



M 



(6) 



where the "time step" r — (3/M. Usually the hamilto- 
nian H is a sum of a kinetic part and a potential one: 
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■H = /C + V. For very small values of r one can use the 
so-called primitive approximation 15| : 



So, for very large values of Af , one has: 
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The density matrix of a system of particles in the 
primitive approximation is given by the following path 
integral: 

p(Ro,Rm;/3) = (AttXt)''^''/^ J ... J dRi...dRM-i 
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-(y(R,_i) + mO) 



The system's partition function writes: 



JJ dRodRA//9(Ro,RAf;/3) 



(9) 



(10) 



Note that this partition function is similar to the parti- 
tion function of a classical system of N polymers. Thus, 
evaluating the average value of any observable of the ini- 
tial many-body quantum system is equivalent to evaluat- 
ing the average value of the observable in the classical 
system described by the partion function Z given by eqs. 
(9), (10). 

So far we dealt with distinguishable particle wave func- 
tions. As we know the bosonic wave functions are symet- 
ric and the fermionic ones are antisymetric. This writes: 



(11) 



where the indexes B and F stand for bosonic / fermionic 
and V is one of the A'^! permutations between the parti- 
cle labels of the N particle the many-body coordinate 
TZ. The sign + corresponds to bosonic systems and 
the sign — to the fermionic ones. One can obtain the 
(anti)symetrization by applying the (anti)symetrization 
operators on the distinguishable particle wave-function: 

$b/f(R) = -j= Y.{±r'^D{V^) (12) 



SO, 



PB/i=^(R,R';/3) = ^^(±)^p,,(R,PR';/3) (13) 



and, replacing pD with the corresponding path integral: 

Pb/f(R, R'; /3) - ^ / - / dRi-dR. 



pi3(R, Ri;/3)...pi5(Rm-i,7^R'; 



(14) 



Now, having the bosonic and fermionic partition func- 
tions of the polymer-like system we can readily perform 
Metropolis Monte Carlo simulations in order to estimate 
average values of various observables. The principle is to 
generate a trajectory in the system's configuration space 
in agreement with the detailled ballanced principle. We 
will not insist here on the simulation since such methods 
are extensively explained elsewhere fl3|. However, it is 
worth noticing that when sampling fermions we deal with 
both positive and negative statistical wheights since the 
average value of any observable writes: 



(0) = 



Ep(-l)^ / dRdR' (R|0|7^RO p(R, PR'; jS) 
Ep(-l)^/dRdR'p(R,7'R';/3) 



(15) 

Therefore, we have the ratio between two differences. 
The problem arises at small temperatures where both 
differences are close to zero so that the statistical fluctu- 
ations increase dramatically. This effect is known as the 
fermion sign problem. However, since we perform simula- 
tions for highly excited nuclear-like systems the fermion 
sign problem shouldn't play here a big role. 

Using the above exposed theoretical ingredients we 
perform Monte Carlo simulations for LJ bosonic and 
fermionic systems. The Lennard- Jones 6-12 (LJ) poten- 
tial writes: 



vo{r) = 4e 



(- 

V r 



(16) 



We use the truncated and long-range corrected version of 
the above potential: v{r) = wo(f) when r < r^] v{r) — 
when r > Tc, corrections being subsequently included in 
order to account for the effect of the neglected tail [T^ ■ 
We took re = L/2, L being the size of the recipient, 
a = 2.55 fm and e — 25.03 MeV. These values are rather 
arbitrary but roughly describe the binding energies of 
various size nuclei. We simulated 16 quantum particles 
placed in a cubic recipient with periodic boundary condi- 
tions. Then, we construct the phase diagram of this sys- 
tem in both bosonic and fermionic quantics. To this aim, 
we estimate pressure versus volume curves at constant 
temperature. Pressure can be easily evaluated starting 
from its canonical definition: P = Td\nZ{/3,V)/dV, 
where Z{/3, V) is the system's canonical partition func- 
tion as defined above. One gets the virial expression for 
pressure: 



P = 



V 3V\^ 



dra 



(17) 



where v{rij) is the particle-particle interaction and () has 
the meaning of canonical average. 

Now let us discuss the results of the simulation. In 
Fig. 1 we represented pressure versus volume curves 
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FIG. 1: Isothermal bosonic and fermionic presure versus vol- 
ume curves corresponding to various temperatures for a 16 
particle system. 



corresponding to both fermionic and bosonic 16 parti- 
cle systems at temperatures ranging from 6 to 10 MeV. 
One can observe the backbendings of these curves, sig- 
natures of a first order liquid-gas phase transition. For 
a given temperature one can observe that all fermionic 
curves are above the bosonic ones. This can be better 
observed in Fig. 2 where bosonic and fermionic pressure 
versus volume curves corresponding to a temperature of 7 
MeV are represented. One further performs Maxwell con- 
structions on all collected pressure versus-volume curve 
and identify the liquid-gas coexistence region for both 
fermionic and bosonic systems. This result is illustrated 
in Fig. 3. There one can observe the phase diagrams 
of both fermionic and bosonic systems. Some interest- 
ing aspects are to be noticed: The fermionic coexistance 
region is wider than the bosonic one. The distance be- 
tween the fermionic and the bosonic curves is larger on 
the denser branch of the phase diagram as expected since 
the quantm interaction between particles is supposed to 
be larger in that region. Finally, the critical point of both 
bosonic and fermionic systems appears to be the same, 
having the value of 10 MeV. This last point is particu- 
larly interesting showing that after the critical point the 
quantum effects are no longer important, after that re- 
gion both systems having a gas-like behaviour. Taking 
into account the small size of the systems (A=16), the 
obtained criticall point temperature, is quite good. In 
infinite nuclear systems the critical temperature is sup- 
posed to be higher, critical temperature being smaller as 
the size of the system decreases 0] . Indeed, from Ref. 
one can deduce a ratio of about 1.58 between the critical 
temperature of an infinite system and the critical tem- 
perature of an A=16 system. One can therefore multiply 
the critical temperature found herein of 10 MeV with the 
factor 1.58 and deduce the critical temperature of an in- 
finite system, T = 15.8 MeV which is a quite realistic 
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FIG. 2: Isothermal bosonic and fermionic presure versus vol- 
ume curves corresponding to a T—7 MeV temperature for a 
16 particle system. 
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FIG. 3: Phase diagram of a bosonic and a fermionic 16 particle 
in a box with PBC conditions system. 



value. 

One can conclude that the Path Integral Monte Carlo 
is a powerfull tool for investigating highly excited nu- 
clear systems. Though the systems under consideration 
are rather small, still valuable information is obtained 
concerning the quantum effects on the systems' phase 
diagram. It was shown that both bosonic and fermionic 
have a common critical point. Moreover, the critical tem- 
perature value of 10 MeV obtainded for an A=16 system 
appears to be realistic, leading via extrapolation to a 
value of 15.8 MeV for infinite systems. 
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